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Abstract 

An essential role in the asteroidal dynamics is played by the mean mo- 
tion resonances. Two-body planet-asteroid resonances are widely known, 
due to the Kirkwood gaps. Besides, so-called three-body mean mo- 
tion resonances exist, in which an asteroid and two planets participate. 
Identification of asteroids in three-body (namely, Jupiter-Saturn-asteroid) 
resonances was initially accomplished by D.Nesvorny and A. Morbidelli 
(1998), who, by means of visual analysis of the time behaviour of res- 
onant arguments, found 255 asteroids to reside in such resonances. We 
develop specialized algorithms and software for massive automatic identi- 
fication of asteroids in the three-body, as well as two-body, resonances of 
arbitrary order, by means of automatic analysis of the time behaviour of 
resonant arguments. In the computation of orbits, all essential perturba- 
tions are taken into account. We integrate the asteroidal orbits on the time 
interval of 100000 yr and identify main-belt asteroids in the three-body 
Jupiter-Saturn-asteroid resonances up to the 6th order inclusive, and in 
the two-body Jupiter-asteroid resonances up to the 9th order inclusive, in 
the set of ~ 250000 objects from the "Asteroids - Dynamic Site" (AstDyS) 
database. The percentages of resonant objects, including extrapolations 
for higher-order resonances, are determined. In particular, the observed 
fraction of pure-resonant asteroids (those exhibiting resonant libration on 
the whole interval of integration) in the three-body resonances up to the 
6th order inclusive is « 0.9% of the whole set; and, using a higher-order 
extrapolation, the actual total fraction of pure-resonant asteroids in the 
three-body resonances of all orders is estimated as w 1.1% of the whole 
set. 

Keywords: Asteroids; Asteroids, dynamics; Celestial mechanics; Resonances 

1 Introduction 

A substantial role of resonances in the dynamics of asteroids became evident 
with the discovery of resonant "gaps" in the asteroid belt by D. Kirkwood in 
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1867. The deepest minima in the distribution of asteroids in the semimajor axes 
of their orbits correspond to the mean motion resonances 2/1, 3/1, 4/1, 5/2, 
and 7/3 with Jupiter. Mean motion resonance represents a commensurability 
between the mean frequencies of the orbital motions of an asteroid and a planet. 
Apart from the mean motion resonances, so-called secular resonances (Murray 
and Dermott, 1999; Morbidelli, 2002), representing commensurabilities between 
the precession rates of the orbits of an asteroid and a planet, are important in 
forming the dynamical structure of the asteroid belt. 

There are two important classes of the mean motion resonances: apart from 
the usual (two-body) mean motion resonances of an asteroid and a planet, an 
appreciable role in the asteroidal dynamics is played by so-called three-body 
mean motion resonances (Murray et al., 1998; Nesvorny and Morbidelli, 1998, 
1999; Morbidelli, 2002). In the latter case, the resonance represents a commen- 
surability between the mean frequencies of the orbital motions of an asteroid 
and two planets (e.g., Jupiter and Saturn): 

mjAj + msAs + m\ « 0, (1) 

where Aj, As, A are the time derivatives of the mean longitudes of Jupiter, 
Saturn, and asteroid, respectively, and mj, mg, m are integers. 

In view of the "overdensity" of the three-body resonances in the phase space 
of the asteroidal motion, Nesvorny and Morbidelli (1998) asserted that "the 
three-body mean motion resonances seem to be the main actors structuring the 
dynamics in the main asteroid belt" . 

Chaotic behaviour, which is often present in the dynamics of celestial bodies, 
is usually due to interaction of resonances (as in any Hamiltonian system, see 
Chirikov 1979), but not always it is known which are the interacting resonances 
that give rise to chaos. It is especially difficult to identify three-body resonances. 
How to distinguish between resonant and non-resonant motions? To solve this 
problem, a "resonant argument" (synonymously "resonant phase" or "critical 
argument") is introduced. It is a linear combination of some angular variables 
of a system under consideration; in the planar asteroidal problem it is given by 

a Pj,Ps,P = m .i-^.i + m s^s +m\+ pjwj + p s zu s + pw, (2) 

where Aj, As, A, rajj, vjs, are the mean longitudes and longitudes of perihelia 
of Jupiter, Saturn, and an asteroid, respectively, and raj, ms, m, pj, ps, p are 
integers satisfying the D'Alembert rule (Morbidelli, 2002): 

mj + m s + m + pj + p s + p — 0. (3) 

If resonant argument (2) librates (similarly to librations of a pendulum), 
the system is in resonance; if it circulates, the system is out of resonance. The 
motion of the system at the border between librations and rotations corresponds 
to the separatrix. Thus the pendulum dynamics provides a graphical model of 
resonance. In a certain sense this model of resonance is "universal" (Chirikov, 
1979). In particular, the motion in three-body resonances can be described in 
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the perturbed pendulum model (Murray et al., 1998; Nesvorny and Morbidclli, 
1998, 1999; Shcvchcnko, 2007). 

An important parameter of a mean motion resonance is its order q, equal to 
the absolute value of the algebraic sum of the coefficients at the mean longitudes 
in the resonant argument: 

q = | to j + TOs + to|. (4) 

The resonant order q is important, because it is the power in which the eccen- 
tricity is raised in the coefficient of the leading resonant term in the expansion 
of the perturbing function (Nesvorny and Morbidelli, 1998). The corresponding 
subresonance width (characterizing also its "strength") is proportional to the 
square root of this coefficient. Thus the value of q determines this important 
property of the leading subresonance. 

In the case of two-body resonances, the role of the resonant order q (defined 
below in Section 6) is analogous: the coefficient of the leading resonant term is 
proportional to e q (Nesvorny and Morbidclli, 1998), where e is the asteroidal 
eccentricity. 

However note that, when there is no strong overlapping of subresonances, 
the resonant order q is not related to the width of the whole resonant multiplet, 
because the separation of subresonances depends solely on the secular preces- 
sion rates of the pericentres (Nesvorny and Morbidelli, 1999); thus the degree 
of overlap (and hence, chaos) in the multiplets is expected to asymptotically 
decrease with the resonant order (Nesvorny and Morbidclli, 1999; Morbidclli 
and Guzzo, 1997). 

One may expect that, generally, broader the leading subresonance of a mean 
motion resonance, greater is the number of objects residing in this mean motion 
resonance. However, no strict correlation exists, due to a competition of various 
dynamical and physical processes, populating or depopulating the resonances. 
We shall discuss this further in more detail. 

In our procedure of resonance identification, described in detail below, we 
limit the set of possible combinations of the integers TOj, to s , to by adopting 
the following conditions: 

q < <7max, (5) 



|toj|, |to s |, |to| < M max , (6) 

where <7, nax = 6 and M max = 8. 

Nesvorny and Morbidelli (1998) used condition (5) (presumably with q max = 
10, as follows from data in table 3 in Nesvorny and Morbidelli 1998). In 
(Nesvorny and Morbidelli, 1999), instead of (6), the following truncation condi- 
tion was used: 

\mj\ + |m s | + |to| < Q 

max (7) 
(see eqs. (29) and (30) and comments on them in Nesvorny and Morbidelli 1999). 
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We identify the three-body resonances in the current motion of asteroids 
with known orbital elements. The limitations of our study are as follows: solely 
the asteroids in the main belt are considered (i.e., the semimajor axes are in the 
range from 2 to 4 AU); solely the three-body resonances with Jupiter and Saturn 
are taken into account; the resonances are considered in the planar problem, i.e., 
the longitudes of nodes in the expression for the resonant argument are ignored; 
the maximum considered order q max of the three-body resonances is set equal 
to 6. 

Our project is intended for the resonance analysis of the orbital data pre- 
sented at the "Asteroids - Dynamic Site" (AstDyS) maintained by A. Milani, 
Z.Knczevic and their coworkers (http://hamilton.dm.unipi.it/cgi-bin/astdys/). 
We take the orbital data for the analysis from this database. Thus the total set 
under analysis contains w 250000 objects. 

The basic purpose of our work is to identify the current three-body reso- 
nances that all the asteroids from the given set are currently involved in. More 
specifically, each object from the set should be put in correspondence to a three- 
body resonance (or none, if there is no resonance). The first attempt of massive 
identification of asteroids in three-body resonances was made by Nesvorny and 
Morbidelli (1998): 255 objects were identified to be in three-body resonances. 
The libration/circulation of the resonant argument for asteroids suspected to 
reside in the resonances was analyzed visually. In our case the data set is much 
greater, and therefore the procedure ought to be completely automatic. Be- 
sides, here we apply a unified bound on the order. This allows one to construct 
a homogeneous identification list for a further statistical analysis. 

To form a general statistical view of the resonant structure of the main belt, 
we also accomplish a massive identification of asteroids in two-body resonances 
with Jupiter, and compare the abundances of asteroids in three-body and two- 
body resonances. 

2 The identification matrix 

As a first stage of the identification process we build an "identification matrix" . 
It consists of two main columns. The first one contains designations of reso- 
nances, and the second one contains the corresponding resonant values of the 
semimajor axis. 

The designations of resonances are given in the notation mjmsm(q) , where 
mj, ms, m are the integer coefficients in the resonant argument (2), and q is the 
resonant order, as defined above. The values of toj, mg, m are given with their 
signs. Thus, examples of this notation look like as follows: 5-2-2(1), 2+2-1(3). 

We construct a set of the resonant arguments for all possible three-body 
resonances up to a fixed order q max in the following way. 

We fix the maximum absolute value M max of each integer mj , m§ , m to be 
equal to g max + 2, where q max is the maximum resonant order. It is assumed 
that 
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mj > 0, gcd(mj,m s ,m) = 1, (8) 

where "gcd" stays for the greatest common divisor. It is set to be equal to 1 
to avoid the higher order harmonics with greater multiplicity. The multiplicity 
is defined as equal to gcd (mj, mg, m). The harmonics with multiplicity greater 
than I are not discernible in our identification procedure, because their argu- 
ments librate simultaneously, though with different amplitudes. (Consider, e.g., 
such resonances as 4 — 2 — 2 and 8 — 4 — 4. The second one, which has multiplicity 
equal to 2, in our procedure is set to be equivalent to the first one.) 

Then we search through all possible combinations of toj , mg , m and identify 
those satisfying the D'Alembert rule (3) and our technical restrictions (5) and 
(6). 

Let us demonstrate how the resonant value of the semimajor axis is calcu- 
lated. According to the definition of the three-body resonance (I), the time 
derivative <jp J;PSiP should be equal to zero. Let us, following Murray et al. 
(1998), assume that w « 0, in the first approximation. Then, for the resonant 
value of mean motion, one has 

n ros = - — (mjAj + m s A s + pjmj + p s m s ). (9) 
m 

Using Kepler's third law, one obtains for the resonant semimajor axis 

(10) 



k N 2 ' 3 



res 



where n les is given by formula (9), and k is the Gauss constant. 

Holman and Murray (1996) and Murray et al. (1998) obtained an approxi- 
mate formula for the precession rate of an asteroid's orbit: 

fi a \ 2 



, , , e*n 3 , (11) 

where \i w 1/1047 is the mass of Jupiter in units of the mass of the Sun, nj 
and aj are Jupiter's mean motion and the semimajor axis of Jupiter's orbit, 
respectively, and 

e =*Z£. (12) 
a.J 

For a, we substitute here the resonant value of the semimajor axis as given by 
Eq. (10). Thus the value of w is calculated. Iterating, one obtains an adequate 
value of a rcs . The second primary column of the identification matrix is filled 
with the resonant values of the semimajor axes, calculated as described, with 
the accuracy of no less than 10~ 3 AU. This accuracy far exceeds the necessary 
one, because we check the asteroids for belonging to a given resonance in a far 
greater neighborhood (±10~ 2 AU) of the computed resonant value of a rcs . 

Nesvorny and Morbidelli (1999) calculated a rcs of the leading subresonances 
(i.e., of the multiplet components 0o,o,-mj-m s -m) °f 19 three-body resonances 
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with the accuracy of 10~ 4 AU: they equated the time derivative of Eq. (2) to 
zero, and, using the values of Aj and As as given by Bretagnon (1990) and w as 
found using the code by Milani and Knezevic (1994), calculated A. All values of 
a ros given in (Nesvorny and Morbidclli, 1999, table 1) agree quite closely with 
the corresponding values of a res calculated here iteratively, as described above, 
for our matrix. The agreement is illustrated in Table 1, where an extract from 
the identification matrix is presented. 



Table 1: An extract from the identification matrix 



mj 


m s 


m 


q 


a rcs (AU) 
(this study) 


a rcs (AU) 
(Nesvorny and Morbidclli, 1999) 


2 


3 


-1 


4 


2.3912 




4 


-2 


-1 


1 


2.3978 


2.3977 


2 


2 


-1 


3 


2.6148 


2.6155 


4 


-3 


-1 





2.6232 


2.6229 


1 


4 


-1 


4 


2.7432 




3 


-1 


-1 


1 


2.7527 


2.7525 


1 


3 


-1 


4 


3.0673 




3 


-2 


-1 





3.0798 


3.0790 


5 


-7 


-1 


3 


3.0925 




5 


-2 


-2 


1 


3.1744 


3.1751 



3 Dynamical identification 

We use the following procedure of dynamical identification. 

First of all, each asteroid's orbit from the adopted set of 249567 objects 
is computed for 10 5 yr. The perturbations from all planets (from Mercury to 
Neptune) and Pluto are taken into account. The hybrid integrator of mer- 
cury6 package (Chambers, 1999) is used. In some cases the orbit9 integrator 
(http://adams.dm.unipi.it/~orbmaint/orbfit/) is used as well, to verify the re- 
sults. The computed trajectories are kept in files for further usage. The trajec- 
tories are output with the time step of 1 yr. 

After integration is over, the objects are taken from the adopted set, and the 
mean value of the semimajor axis is computed for each object. Using this value, 
a set of preliminary resonant arguments cr rcs is found in the identification matrix. 
Each argument cr rcs is then analyzed on the presence of libration/circulation, 
using the computed trajectory of the object. 

We distinguish two types of resonant libration: pure and transient. By 
definition, the libration is pure, if it lasts during the whole time interval of 
integration, i.e., 10 5 yr. An example of such libration is given in Fig. 1, where the 
time behaviour of the resonant argument, alongside with the orbital elements, 
is demonstrated for asteroid 463 Lola, resonance 4 — 2 — 1. 
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The libration is defined as transient, if circulation appears at any time during 
this interval. An example is given in Fig. 2; this is the case of asteroid 490 Ver- 
itas, resonance 5 — 2 — 2. In Fig. 2, it may seem unusual that Veritas exhibits 
circulation of the resonant argument while the semimajor axis remains almost 
constant (especially in comparison with Fig. 9, which is considered below and 
which shows the orbital elements and resonant argument of 1915 Quetzalcoatl, 
residing in the two-body transient resonance 3/1). This difference is explained 
by a large difference in the widths of the resonances: e.g., at e = 0.1 the width 
of resonance 3/1 is ~ 5 times greater than that of resonance 5 — 2 — 2 (see fig. 1 
in Morbidelli and Nesvorny 1999). Moreover, 1915 Quetzalcoatl has a very large 
eccentricity (~ 0.6-0.8). At e = 0.1, the half-width of the 5 — 2 — 2 resonance is 
w 0.002 AU (sec table 1 in Nesvorny and Morbidelli 1999 and fig. 1 in Morbidelli 
and Nesvorny 1999). For Veritas, the eccentricity is ~ 0.06 (see Fig. 2), hence 
the half- width is even smaller. This makes the shift in a, when a is circulating, 
almost imperceptible, especially in the digitally unfiltcred a. When a is filtered, 
such shifts look more obvious; see fig. 2 in (Knezevic et al., 2002) and fig. 3 
in (Tsiganis et al., 2007). Note that, when the circulations are short-term, the 
shifts in a can be imperceptible even in the digitally filtered element: see fig. 1 
in (Nesvorny and Morbidelli, 1998), where the time behaviour of Veritas in the 
digitally filtered elements is shown. 

To distinguish between transient-resonant and non-resonant behaviours we 
introduce a technical parameter: the resonance minimum time, which we set 
to be equal to 20000 yr. Roughly speaking, this is the minimum time interval 
of resonant librations, for an asteroid to be regarded as resonant. An exact 
definition is given below. 

Nesvorny and Morbidelli (1998) identified resonant behaviour visually and 
used more subjective criteria, according to which the asteroid is resonant if 
"(1) the corresponding resonant angle shows evident librations during the in- 
tegration time span or (2) the resonant angle circulates with a period longer 
than several thousand years" (Nesvorny and Morbidelli, 1998). In fact the sec- 
ond criterion threshold is analogous to our resonance minimum time, because 
our automatic procedure measures the time intervals of libration as such when 
no circulation is present, as described in detail below. Our criterion thresh- 
old, being larger, is more restrictive; however, in both cases the threshold is 
much greater than the timescales of circulation far from resonance; the latter 
timescales are of the order of asteroidal orbital periods. 

The complete procedure of identification is as following. When computing 
each trajectory, the resonant argument value for the guiding subresonance is 
calculated at each step in time, and these values are written in a file, in function 
of time. After the trajectory computation is over, the time behaviour of the 
resonant argument is analyzed. If a current value of the resonant argument is 
different from its previous value by more than 2n 7 this means that there is a 
break, which takes place if there is either circulation or apocentric libration. To 
distinguish between these two cases, the same procedure is repeated, but with an 
artificial shift of the resonant argument. This shift is equal to ir. Then, if there 
is apocentric libration, it turns into pericentric one, and there are no breaks. 
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At the beginning of the trajectory analysis, two variables are initiated: the first 
one is the duration of current libration, and the second one is the total time of 
libration. If circulation whenever starts, the first variable value is added to the 
second one and the first variable is reset to zero. At the end of the analysis, it is 
checked whether the second variable (total time of libration) is equal to the full 
time of computation (10 5 yr). If yes, then the asteroid is regarded to be in pure 
resonance. If no, but the second variable value exceeds the resonance minimum 
time, then the asteroid is in transient resonance. If the second variable value 
is less than the resonance minimum time, then the asteroid is regarded to be 
non-resonant. The adopted pure/transient division of resonances is thus mostly 
technical, because it depends on the chosen time interval of computation. 

We use solely the direct method of identification, i.e., the resonant argument 
is checked on the subject of libration. Any secondary and/or auxiliary criteria, 
such as the semimajor axis behaviour, the location of objects in the proper 
"semimajor axis - eccentricity" plane, the Lyapunov exponents are left for future 
explorations, — with a partial exception in the latter case, see Section 5. 

In our identification procedure, it is formally possible that an asteroid might 
be identified as belonging to two or even more resonances, because all resonances 
in a rather broad neighbourhood in semimajor axis (±I0~ 2 AU, as mentioned 
above) in the identification matrix are scanned on the subject of libration. If two 
three-body resonances are close enough in semimajor axis, and especially if they 
also overlap with a strong two-body resonance, an asteroid may intermittently 
diffuse from one to another resonance; a small number of such objects have 
been identified, as belonging to two resonances. In our lists, each object of this 
kind has been finally attributed to the resonance where it stayed for a longer 
time, to avoid complication of statistics. Analysis of such objects (we call them 
"rogue-resonant" asteroids) will be given elsewhere. 

Besides, since the data used is unfiltered, some misidentifications may take 
place when the libration amplitude is high. To estimate the formal statistical 
error of our procedure, we randomly chose 300 asteroids identified as pure- 
resonant and visually checked the libration of their resonant arguments. It 
turned out that only 4 objects (1.3% of the set) were misidentificd. Thus the 
statistical error of the procedure is « 1%. 

All identified pure three-body-resonant asteroids, grouped according to as- 
sociation to a given resonance, are listed in the Appendix A. The top ten res- 
onances, that are most "populated" , are listed in Table 2. The last column of 
Table 2 contains analytical estimates of the resonance width (at the asteroidal 
eccentricity e = 0.1), according to (Nesvorny and Morbidelli, 1999, table 1). 
One can see that, generally, broader the resonance, greater is the number of ob- 
jects residing in it. However, this tendency is not strict. The reason is that the 
dynamics here is strongly interrelated with physics: e.g., a collisional disinte- 
gration of an asteroid can strongly increase abundance of objects in a particular 
resonance, thus disturbing the expected correlation. 

For each resonance, statistics on the asteroids in pure and transient librations 
have been calculated. The statistical results are summarized in Table 3. The 
fraction of asteroids in three-body resonances (transient plus pure) turns out to 
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Figure 1: The orbital elements and resonant argument of 463 Lola. Resonance 



be ~ 4.4% of the whole set. This is rather close to the value 4.6% (« 1500 in 
the set of « 32400) found by Nesvorny and Morbidclli (1998) to serve as a lower 
bound for the relative number of resonant asteroids. As follows from Table 3, 
the fraction of asteroids in pure three-body resonances turns out to be « 0.94% 
of the whole set. 

The third column of Table 3 contains a prediction for the numbers of as- 
teroids residing in resonances of all orders; this subject is discussed in the next 
Section. 

Let us consider the location of three-body-resonant asteroids in the "semi- 
major axis — eccentricity" plane. For each asteroid in a given resonance we 
calculate the average values (over the whole time interval of integration) of the 
semimajor axis and eccentricity and plot these values in the "a-e" plane. We 
have accomplished this procedure for each resonance with known structure of 
the resonant multiplet, as calculated by Nesvorny and Morbidelli (1999), so that 
the separatrix of the leading subresonance could be drawn. An example of such 
a plot is given in Fig. 3. It is clear that the percentage of "outliers" (objects 
out of the separatrix cell) is zero. This confirms the good accuracy of the ac- 
complished identification procedure in the case of this particular resonance. A 
detailed study of the "a-e" plots will be given elsewhere. 

We have cross-checked the lists of resonant objects identified in our study and 
in (Nesvorny and Morbidelli, 1998). It turns out that the number of resonant 
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Figure 2: The orbital elements and resonant argument of 490 Veritas. Resonance 
5-2-2. 

objects listed by Ncsvorny and Morbidclli (1998) but not identified as resonant 
in our study docs not exceed 1% of the list of Nesvorny and Morbidclli (1998). 
This confirms that the differences in the methodologies of identification in the 
two studies do not play any significant role in what concerns the reliability of 
results. 

4 Expected abundances of asteroids in high- 
order three-body resonances 

The obtained list of resonant asteroids is obviously not complete, due to the 
limitations of identification criteria. First of all, the resonant order q < 6. 
Presumably, there is a lot of objects in resonances of higher orders. To take 
account of them, let us analyze the asteroid distribution in the resonant order 
q. The necessary data as derived from the results of the identification procedure 
are presented in Table 4. 

The constructed differential distributions (histograms) in q are shown in 
Fig. 4 for the case of transient plus pure resonances, and in Fig. 5 for the case 
of pure resonances. 

As an approximating function, we have chosen the power law 
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a, AU 

Figure 3: The asteroids identified to be in resonance 5 — 2 — 2: location in 
the "a-e" plane. Solid curve: the separatrix of the leading subresonance, as 
calculated in (Nesvorny and Morbidelli, 1999). 

Table 2: The top ten most populated three-body resonances Sun-Jupiter- 



asteroid 


mj m s 


m 


Ores 


Trans. + pure, 


Pure, number 


Res. width (AU) at e = 0.1 






(AU) 


number of objects 


of objects 


(Nesvorny and Morbidelli, 1999) 


5 -2 


-2 


3.1744 


699 


182 


0.0056 


4 -2 


-1 


2.3978 


688 


595 


0.0024 


3 -2 


-1 


3.0798 


621 


134 


0.0045 


3 -1 


-1 


2.7527 


540 


203 


0.0019 


2 2 


-1 


2.6148 


470 


34 


0.00015 


4 -3 


-1 


2.6232 


455 


90 


0.0009 


2 3 


-1 


2.3912 


343 


56 




5 -7 


-1 


3.0925 


314 


48 




1 3 


-1 


3.0673 


300 


52 




1 4 


-1 


2.7432 


284 


47 





N = aq b , (13) 

where a and b are two fitting parameters; b < 0. We have also tried the exponen- 
tial law N oc exp(cg) (where c < 0), but it has turned out to be inappropriate, 
the statistical significance of fitting being very low. 

Since we are interested in the tail behaviour of the distributions, we have 
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Table 3: Asteroids in three-body resonances, statistics 



NM98* This study, This study, prediction 







<?max — 6 


for <7 max = oo 


The whole set of objects 


5400 


249567 


249567 


Objects with integrated orbits 


836 


249567 


249567 


Transient+pure-resonant objects 


255 


11039 




The same, fraction of the studied set 


4.7% 


4.4% 




Pure-resonant objects 




2338 


2854 


The same, fraction of the studied set 




0.94% 


1.1% 



*Nesvorny and Morbidelli (1998). 



Table 4: Abundances of asteroids in three-body resonances in function of q 



q 





1 


2 


3 


4 


5 


6 


Transient+pure 


1371 


2884 


1367 


2256 


1394 


997 


770 


Pure 


356 


1119 


194 


251 


214 


117 


87 



used the data for q > 1, ignoring the specific case q — 0. In the case of all 
resonant (transient-resonant plus pure-resonant) asteroids (see Fig. 4) we find 
a = 2821 ± 454, b = -0.57 ± 0.18; and the correlation coefficient R 2 = 0.71. As 
soon as |6| < 1, the predicted number of asteroids in high-order resonances is 
formally infinite; in practice this means that it can comprise up to ~ 100% of 
the whole set. 

In the case of pure-resonant asteroids (see Fig. 5) we find a = 1095 ± 103, 
b = -1.64 ± 0.28; R 2 = 0.94. As soon as |b| > 1, the predicted total number 
of asteroids in high-order resonances is finite. Using Eq. (13), the number of 
objects with q > 7 in the studied set is estimated to be equal to w 516, or 
« 22.1% of the identified number (2338). Therefore, the predicted total number 
of asteroids in pure three-body resonances is estimated as 2854, constituting 
1.1% of the whole set (249567). 

A note of caution is in order here. One has to admit that the fits made in 
this Section are based on few points, and, therefore, any statistical predictions, 
made with these formulas, are uncertain. What is more, one cannot even expect 
to find smooth distributions and/or strict correlations in this field of research, 
where dynamics is strongly interrelated with physics: the abundances of objects 
in resonant groups are regulated by various processes, e.g., by collisions and 
the Yarkovsky effect. Interactions with asteroidal families are important. Of 
course, higher order resonances should be directly analyzed in the future. 

Another complicating factor is that the subresonance overlap (and hence, 
the degree of chaoticity) in the multiplets is expected to decrease asymptoti- 
cally at q 3> 1. The reason is that the subresonance typical width scales with 
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Figure 4: Circles: the distribution of all resonant (transient-resonant plus pure- 
resonant) asteroids in resonant order q > 1. Solid curve: the power-law fitting. 



q as ~ e q (where e is the eccentricity), whereas the subresonances separation 
(determining the multiplct width) remains basically constant, depending solely 
on the secular precession rates of the pericentres; thus the overlap/interaction 
of subresonances in the multiplets decreases asymptotically at q ^> 1 (Nesvorny 
and Morbidclli, 1999, p. 268-269). When the ratio of the subresonances sepa- 
ration to the width of the leading subresonance is much greater than one, the 
separatrix chaotic layers are exponentially thin with this ratio, see (Chirikov, 
1979; Shcvchcnko, 2008, 2011). It is straightforward to suppose that, if one 
fixes the asteroidal eccentricity, three-body resonances with increasing order 
become basically regular (Nesvorny and Morbidclli, 1999). Thus, if the reso- 
nances are populated uniformly in the eccentricity, the power-law extrapolation 
for the high-order transient-resonant populations fails beyond some order. Only 
if this critical order is high enough, the extrapolation-based estimates might be 
appropriate. 
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Figure 5: Circles: the distribution of pure-resonant asteroids in resonant order 
q > 1. Solid curve: the power-law fitting. 

5 Lyapunov exponents in three-body reso- 
nances 

The AstDyS database provides information on the maximum Lyapunov expo- 
nents for almost all asteroids contained in it. The provided values are computed 
on time intervals of 2 mln yr, i.e., on time intervals 20 times greater than that 
we use in our identification procedure. It is instructive to check how the data 
on Lyapunov exponents correlate with the pure/transient division (adopted in 
our study) of the resonant asteroids. 

In the transient case the motion is expected to be chaotic, and in the pure 
case to be essentially regular, because transitions from libration to circulation 
and vice versa are inevitable in a separatrix chaotic layer of resonance. However, 
the adopted relatively short integration time may imply that many "pure"- 
resonant asteroids become transient on longer timescales. 

First of all, let us build distributions of resonant asteroids in L (the maximum 
Lyapunov exponent). The constructed differential distributions (histograms) 
are presented in a single plot in Fig. 6. The Lyapunov exponents are given 
in units of (mln yr) -1 . (Note that in the AstDyS database they are given in 
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units of yr -1 .) N is the number of objects in the interval (L, L + AL), where 
AL = 10. The grey histogram is for the asteroids in transient resonances, and 
the black one is for the asteroids in pure resonances. The histograms are cut off 
at L = 500, because the objects with greater values of L are rare; they are all 
transients. 

From Fig. 6 it is clear that the transients have a much more extended dis- 
tribution in L, in comparison with the pure-resonant objects, as expected. This 
difference is uniform in all resonant groups, as can be directly seen from Fig. 7, 
where the maximum Lyapunov exponents are plotted versus the semimajor axis 
a: the positions of open dots (representing transients) extend to much greater 
heights, in comparison with the pure-resonant objects, in all resonant groups. 
Note that this plot is built without any cut-off in L, i.e., all objects are shown. 
It is evident that only transients have very large Lyapunov exponents, corre- 
sponding to Lyapunov times as small as w 570 yr. 

10000 
1000 
N 100 
10 

1 

L 



Figure 6: The distributions of resonant asteroids in the value of maximum 
Lyapunov exponent. Grey: asteroids in transient resonances. Black: asteroids 
in pure resonances. 

Among all identified pure- resonant objects, the average L is 34.3 (mln yr) -1 , 
and among all identified transients it is 49.7 (mln yr) -1 . The Lyapunov times 
are s=s 29200 and s=s 20100 yr, respectively. Among all objects of the AstDyS 
database with measured Lyapunov exponents, the average L = 20.8 (mln yr) -1 , 
corresponding to the Lyapunov time w 48100 yr. Thus both pure-resonant and 
transient-resonant objects turn out to be more chaotic than a typical asteroid. 
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Figure 7: The maximum Lyapunov exponents versus the semimajor axis a. 
Open dots: asteroids in transient resonances. Black dots: asteroids in pure 
resonances. 



The fact that many pure-resonant asteroids have definitely non-zero Lyapunov 
exponents signifies that they become transient on the timescales longer than the 
adopted integration time (10 5 yr). This underlines the conditional character of 
the "pure/transient" technical classification adopted in our study. 

6 Two-body resonances with Jupiter 

To form a more general statistical view of the actual resonant structure of the 
main belt, it is instructive to compare the abundances of asteroids in three-body 
resonances with the abundances of asteroids in two-body resonances. For this 
purpose, we have performed a procedure of identification of asteroids in two- 
body resonances with Jupiter, analogous to that described above for the case of 
three-body resonances. In this identification procedure, we assume circular and 
zero inclination orbits of perturbing planets. 

Taking into account the D'Alembert rule, the resonant argument for the 
resonance of order q is defined by the following formula (Murray and Dermott, 
1999; Morbidelli, 2002; Gallardo, 2006): 
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a = (p + q)Xj — pX - qw, 



(14) 



where Aj and A are the mean longitudes of Jupiter and an asteroid, respectively, 
and w is the longitude of perihelion of the asteroid; q is the resonant order, p 
is integer. 

Of course, this setting of the identification problem is a rather simplified one. 
When real resonant asteroids are considered, the "circular" subresonance term 
is not expected to dominate universally; what is more, the resonant asteroid 
might reside in another component of a multiplet (see Holman and Murray 
1996; Murray and Holman 1997). We regard the adopted approach as a first 
approximation, and leave considerations of all possible subresonances for the 
future. 

The resonant value of the semimajor axis of an asteroidal orbit is given by 



where aj is the semimajor axis of Jupiter's orbit, and ^ is the mass of Jupiter 
in units of the mass of the Sun (see, e.g., Murray and Dermott 1999; Gallardo 



The identification matrix is constructed in the way analogous to that de- 
scribed in Section 2. Similar to the case of three-body resonances, it is assumed 
that gcd(p, q) = 1. Besides, 1 < p < 11. 

The resonant order range is limited to < q < 9. Recall that in the case of 
three-body resonances we have set < q < 6. Choosing the upper bounds with 
the difference equal to 3 allows one to adequately compare identification statis- 
tics for two-body and three-body resonances, because "a three-body resonance 
of a given order q should have roughly the same strength as a usual resonance of 
order q + 3 for eccentricity of about 0.05-0.10" (Nesvorny and Morbidelli, 1998). 
Such a trick was used in Nesvorny and Morbidelli (1998) for similar purposes. 

Figs. 8 and 9 show the time behaviour of the resonant argument and or- 
bital elements of 1915 Quetzalcoatl and 190 Ismene, as two typical examples. 
The top ten most "populated" resonances, identified in our study, are listed in 
Table 5. All identified pure two- body- resonant asteroids, grouped according to 
association to a given resonance, are listed in the Appendix B. 1 

The resulting statistics of identified objects arc listed in Tables 6 and 7. As 
follows from Table 6, a half of all identified asteroids in pure two-body resonances 
are Trojans (1669/3132 « 53%). Pure Trojans plus pure Hildas constitute 
« 85% of all asteroids in pure two-body resonances in our set. In Table 7, the q 
dependence of the resonant asteroid abundances is presented. The dependence 
is obviously irregular and does not permit any smooth decay approximation. 
Especially one should point out the negligible asteroid abundances at q = 3 and 



1 One may wonder why 279 Thule is absent in the 4/3 entry. The matter is that our 
automatic procedure identifies it as being in transient resonance, not in pure one, because the 
resonant argument sometimes goes out of the range (— n, +tt). 




(15) 



2006). 



q = 9. 
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Figure 8: The orbital elements and resonant argument of 190 Ismene. Pure 
resonance 3/2. 



Table 5: The top ten most populated two-body resonances with Jupiter. The 
main belt and Trojans 



p + q 


P 


q 


Ores 

(AU) 


Transient 
+ pure 


Pure 


Notes 


1 


1 





5.2043 


1670 


1669 


Trojan swarms 


3 


2 


l 


3.9716 


1033 


1007 


Hilda group 


11 


5 


6 


3.0766 


351 


49 




2 


1 


1 


3.2785 


274 


227 


Griqua family 


8 


3 


5 


2.7063 


232 


20 




7 


2 


5 


2.2576 


208 


20 




10 


3 


7 


2.3322 


165 


15 




11 


4 


7 


2.6514 


165 







9 


4 


5 


3.0309 


163 


55 




3 


1 


2 


2.5020 


63 


32 


Alinda family 


7 


3 


4 


2.9583 


51 
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It turns out that in transient plus pure resonances the asteroids are ~ 2.5 
times more abundant in three-body resonances than in two-body resonances 
(11039 versus 4450); and in pure resonances the abundances are in the ratio 
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Figure 9: The orbital elements and resonant argument of 1915 Quctzalcoatl. 
Transient resonance 3/1. 



Table 6: Asteroids in two-body resonances with Jupiter, statistics. The main 
belt and Trojans 





Number 


Number 




of resonances 


of objects 


Transient + pure 


21 


4450 (1.78%) 


Pure 


16 


3132 (1.25%) 



Table 7: Abundances of asteroids in two-body resonances in function of q 



q 





1 


2 


3 


4 


5 


6 


7 


8 9 


Transient + pure 


1670 


1310 


63 


15 


57 


603 


352 


342 


37 1 


Pure 


1669 


1236 


32 


5 


19 


95 


49 


23 


4 



« 3 : 4 (2338 versus 3132). However, if one excludes Trojans and Hildas, the 
abundance of three-body-resonant asteroids becomes overwhelming: in the case 
of transient plus pure resonances the ratio of observed abundances of three-body- 
resonant objects and two-body-resonant objects becomes equal to 11039/1747 w 
6.3; and in the case of pure resonances the ratio becomes equal to 2338/456 w 
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5.1. 

If Trojans and Hildas arc not excluded, the fraction of pure-resonant aster- 
oids among the transient+pure-resonant asteroids is w 3.3 times greater in the 
two-body case than in the three-body one (« 70.4% versus « 21.2%); thus it 
might seem that the two-body resonances are much more "pure" on average. 
However, if one excludes Trojans and Hildas, the fractions become comparable: 
26.1% versus 21.2%. 

One should outline that the comparative analysis given here is merely sta- 
tistical and, thus, formal. However, such an analysis might provide a necessary 
preliminary stage for a deeper study of the comparative role of two-body and 
three-body resonances; such a study should comprise consideration of the physi- 
cal (collisions and the Yarkovsky effect) and dynamical (transport and diffusion) 
processes, which lead to populating or depopulating the resonances. 

7 Conclusions 

1. We have identified the resonant objects (the objects residing in three-body 
resonances with Jupiter and Saturn in the main asteroid belt) in the set 
of all numbered asteroids in the AstDyS database. This set comprises 
249567 asteroids catalogued up to the date of April, 2011. The list of all 
asteroids identified as residing in pure three-body resonances is given in 
Appendix A. 

2. The fraction of asteroids in three-body resonances (transient plus pure) up 
to the 6th order inclusive turns out to be « 4.4% of the total studied set of 
249567 asteroids. The fraction of asteroids in pure three-body resonances 
of the same orders turns out to be « 0.94% of the total studied set. 

3. The top three most populated three-body resonances are: 5 -2 -2 (contain- 
ing 699 transient+pure-resonant asteroids), 4 -2 -1 (688 transient+pure- 
resonant asteroids), 3 -2 -1 (621 transient+pure-resonant asteroids). For 
the pure-resonant asteroids, the "top three" resonances are: 4-2-1, 3-1 
-1, and 5 -2 -2, containing 595, 203, and 182 objects, respectively. 

4. Using a high-order extrapolation (in the form of a power law) of the q 
dependence of the number of identified resonant objects, the actual total 
fraction of asteroids in pure three-body resonances of all orders is esti- 
mated as « 1.1% of the whole set. In what concerns the case of transient 
three-body resonances, the situation is much less certain, because the 
power-law extrapolation diverges. 

5. We have also identified all objects residing in two-body resonances (of 
order < q < 9) with Jupiter in the main asteroid belt, taking the same 
database of asteroids. The list of all asteroids identified as residing in pure 
two-body resonances is given in Appendix B. 
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6. The half of all identified asteroids in pure two-body resonances are Tro- 
jans (« 53%). The pure Trojans plus pure Hildas constitute « 85% of 
all asteroids in the pure two-body resonances. The q dependence of the 
two-body resonant abundances is clearly irregular and does not permit 
any smooth decay approximation. Especially one should point out the 
negligible asteroidal abundances at q — 3 and q = 9. 

7. In the transient plus pure resonances, the identified asteroids are w 2.5 
times more abundant in the three-body resonances than in the two-body 
resonances; and in the pure resonances the abundances are comparable. 
However, if one excludes Trojans and Hildas, the abundance of three-body- 
resonant asteroids becomes overwhelming. What is more, taking into ac- 
count extrapolated abundances in high-order resonances may substantially 
increase this overwhelming domination. Thus our analysis quantitatively 
verifies the assertion by Nesvorny and Morbidclli (1998) that "the three- 
body mean motion resonances seem to be the main actors structuring the 
dynamics in the main asteroid belt" . 

8. We would like to point out that our results confirm the general concept 
of Molchanov (1968, 1969) on the omnipresence of resonances in the Solar 
system, however at a new level of understanding of this phenomenon. 
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Appendix A 



Asteroids in pure three-body resonances with Jupiter and 
Saturn 

13- 1 (3): 95 423 1102 3585 11456 17433 22499 27767 34493 41164 42441 
46218 47896 48131 59909 62884 63667 72791 73569 74425 76056 78202 82703 
91809 105839 112669 117065 117344 129360 138837 140653 141042 145390 
155877 159321 161779 166668 186951 195633 195930 199688 200344 201223 
201875 209429 209437 209473 210159 218079 226525 232332 232494 238387 
239477 240226 

14- 1 (4): 9835 20596 31058 36595 37004 37611 41056 41640 50350 
51002 55069 57488 58002 62632 66629 67442 71024 82600 83471 84938 87318 
90368 93699 97281 107076 119544 122676 126319 130159 147292 158732 160021 
166433 166516 167947 168639 178642 186262 188442 192877 212621 219368 
226015 232815 235341 241349 249070 

16-1 (6): 10482 12785 14833 16967 18960 25508 31374 33523 34780 
39875 48030 53741 67456 74729 79395 82094 95170 105575 105721 130930 
147022 180268 194048 199297 203374 210376 226238 230570 232392 233526 
245785 

2 1 -1 (2): 1149 1670 38992 67811 105476 167454 229942 241160 

2 2 -1 (3): 70 194 258 839 923 995 16233 28706 37158 48104 48110 
48127 48418 48532 51196 62082 67329 71407 84437 91337 93241 117209 122690 
126001 133792 154671 173752 180242 186145 212288 224957 232199 238618 
243865 

2 3-1 (4): 7553 10368 39930 42942 48007 48028 48039 48105 48111 
48118 48120 48121 48133 48509 48511 69468 70514 72106 74316 82341 89577 
90015 90354 90546 112351 113728 115283 119372 122007 125642 129963 130692 
150940 153684 157179 159616 161644 162752 167044 180486 180491 180493 
180494 186531 188874 189427 197995 203427 216047 220119 222186 236300 
239888 243110 244652 248126 

2 4-1 (5): 9945 12169 15628 32302 32898 33946 38779 38993 39176 
48014 48033 48103 48107 59305 75243 90851 99504 120401 138965 143630 
168287 180218 211512 217139 218091 239263 

3 -2 -1 (0): 9864 10173 10926 17320 22785 26986 29849 29869 34592 
36671 40940 41403 45312 47919 54497 56883 58847 58914 62989 64773 71551 
73469 73570 73573 73798 73952 75981 76535 78756 82905 83011 83055 83103 
83698 83972 84448 84640 90097 94202 97229 106316 112360 112384 112757 
113085 113337 113402 113904 115445 123903 124012 127061 128255 128304 
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129408 130325 131350 132895 132912 133447 133714 135218 136290 138750 

140372 140655 143617 144889 152112 157195 158147 160008 160150 164493 

166768 166873 167567 169214 169298 170112 172023 173990 175484 182584 

186114 189942 190125 195672 195779 195797 196016 196634 198642 199662 

201338 201368 202406 203309 203815 204167 207764 209436 209446 209457 

209458 209482 209487 209488 210167 210978 214804 216414 216642 217332 

218433 221358 221408 223303 226565 227722 228301 231654 234751 234855 
234957 235444 235958 240046 243227 243296 246676 247045 248205 249330 

3 -1 -1 (1): 564 947 1427 2042 2465 3534 4106 4426 5701 8068 9447 
13085 13357 14169 14368 18302 19901 20118 20694 22146 22562 24011 24032 
24190 25314 25891 27827 31192 33071 34104 35926 37025 41157 41868 42743 
44362 44458 45486 48038 48539 50298 50954 54089 54405 54429 55893 56860 
57161 57879 59296 59681 60160 61974 61980 62518 63423 64937 65543 67013 
67028 67349 69300 69375 71281 71401 71423 71428 74215 74225 74314 76226 
76238 76387 78768 79443 81841 81883 81927 83388 83935 87749 87893 88110 



88173 ! 


39354 90647 93212 93749 95486 96994 


97230 


112181 


114698 


115643 


117401 


117668 


118375 


118835 


120109 


122501 


122720 


123213 


123891 


126305 


127119 


127499 


131131 


135585 


138054 


138096 


139828 


139925 


140586 


140605 


140891 


143738 


146483 


146507 


147367 


147370 


148921 


148949 


148965 


149535 


152057 156621 


157792 


159133 


161375 


162113 


164762 


170012 


170043 


170655 


171450 


172327 


172703 


172908 


174418 


177037 


179311 


180571 


180686 


181134 


181576 


181761 


181872 


182146 


185159 


185934 


188312 


192832 


195353 


196624 


198484 


198567 


198813 


199291 


200275 


201564 


201633 


201847 


202295 


203717 


204124 


204800 


206078 


206622 


206654 


207269 


207302 


207661 


208908 


211245 


212444 


213406 


217444 


217785 


217790 


218768 


218955 


219401 


220218 


220322 


222347 223648 


223656 


223845 


225299 


227701 


229534 


230119 


231697 


233422 



234345 235689 239615 241729 242054 244969 245373 245400 247051 

3 1-1 (3): 1705 6581 10851 15431 22249 36579 41431 44783 48124 
50081 52125 52586 53090 58185 59601 64614 66163 68674 68996 69393 69636 
74897 79142 80390 82322 88315 90730 97992 105977 114325 115219 118117 
119371 125682 129748 129751 131646 134132 137237 139138 145918 153401 
162608 165365 176110 177795 179798 180412 197692 205613 213218 213525 
217478 220043 226224 232412 235665 240796 

3 3-2 (4): 5413 10656 11653 21807 29946 41152 48515 48516 106548 
106707 121883 138820 149695 182619 185797 209431 209440 209464 209469 
209478 209483 209494 223389 

3 4-2 (5): 28984 32932 33778 36169 39340 42104 56894 113850 116063 
167477 168326 197436 209435 209439 209443 225316 

3 4 -1 (6): 209454 

4 -4 -1 (1): 52667 107282 129662 157240 163464 198658 
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4 -3 -1 (0): 3870 7456 7530 19700 21151 22085 22686 22904 23678 
23960 33719 35565 37541 37927 39547 39910 42708 48130 48535 50842 57528 
61266 61284 61693 66867 68463 73206 75281 76150 76180 78509 81286 90130 
90690 90721 93932 94135 96877 96966 99102 105019 105230 112405 113756 
113933 119525 126653 126750 126821 132369 135793 138256 139761 144051 
148896 149892 150968 158721 162842 165625 168566 169085 171254 172708 
178376 179009 182622 182698 186902 187380 195541 196456 205418 210213 
212226 214528 215264 221172 222747 227920 229825 236640 238272 240937 
241763 241766 241991 242342 246382 248308 

4 -2 -1 (1): 463 2348 2487 2640 2791 3048 3293 3699 3716 3733 3865 
4684 5156 5323 5747 6167 6193 6947 7133 7831 8044 8166 8772 9777 10009 



10058 


10092 


10218 


10371 


10701 


11011 


11256 


12065 


12097 


12161 


12177 


12299 


12510 


12674 


12730 


12985 


13002 


13024 


13127 


13159 


13637 
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244214 244752 244778 246097 246158 246168 247378 247851 248673 

4 -1 -1 (2): 443 2175 2440 10364 20912 29435 30321 30916 53011 65455 67555 
67699 76876 98342 98448 106572 116738 121837 124515 127569 129739 137412 
141493 141594 143599 147582 149233 151938 171702 177786 178533 203286 
210206 214109 214530 215602 217517 219155 224790 231602 242502 

4 1 -2 (3): 116793 

4 1 -1 (4): 5066 28441 141359 166908 209455 

4 3-2 (5): 37957 62712 81077 81561 115442 126605 170321 180727 
209153 213348 226790 

5 -7 -1 (3): 19282 22198 48500 58123 58525 59375 82464 82946 83299 
84232 84696 87576 97536 105807 113343 118397 121291 130281 138591 141285 
149124 155580 161209 161260 161771 168264 172761 179500 181007 189914 
195767 195867 197461 201308 206843 208354 208477 209442 209444 209463 
209471 209495 212905 217259 224768 237985 238056 242698 

5-7 5 (3): 23539 26521 30985 71319 73466 76454 96142 106325 141402 
145415 231035 240779 

5 -6 -1 (2): 5534 7176 14090 14941 17149 21960 29351 29871 30477 
35153 35575 39299 41108 50657 53321 56310 62238 62810 68099 68723 69897 
70592 70951 97743 105868 114380 121278 122957 123365 125842 132699 133224 
139562 146555 147789 150481 153062 164540 170742 176979 177535 178321 
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195294 201559 221246 224675 230944 237117 243700 249153 
5 -6 5 (4): 43685 219393 

5 -4 -1 (0): 5277 9459 12814 13065 13757 16424 19419 22792 23520 
23657 27100 27804 29650 30328 30679 30959 31439 32804 33412 35718 35991 
36418 36424 38097 39810 42505 43479 44565 45181 45842 45888 46445 47258 
48010 51708 51709 52841 52955 53791 55705 59106 59115 60633 61457 65886 
66194 66269 66295 67385 67829 68668 69573 71753 71909 72327 75046 79865 
80394 80850 88936 89454 92422 92524 97837 111991 112301 114347 115038 
117889 118180 118353 121107 124291 124903 125165 125652 129793 131610 
131808 137413 137578 137760 139104 139153 142094 145952 147019 153446 
155114 155635 158205 160879 161918 163033 165230 166019 170983 171793 
172567 172887 177239 179648 183044 183629 193928 194029 194089 194146 
194297 201491 202894 205754 207455 208634 211549 212425 215432 216611 
217155 217399 221620 222939 226260 229061 230577 232383 236196 238076 
245780 

5 -3 -1 (1): 7882 29998 43935 106173 107454 207395 220455 221522 

5 -2 -2 (1): 490 744 755 786 818 1072 1209 1438 1701 1731 2039 2142 
2164 2250 2492 2587 2666 2731 2757 2863 3460 4152 6830 9026 9163 11973 
12762 12912 14360 16968 17165 17216 18202 21223 21547 21576 22097 25681 
29519 29881 29931 31006 31748 32416 33225 33997 34070 35183 35660 38250 
38507 40718 45587 46095 47652 47949 51926 52257 52693 53351 55682 57870 
59160 62551 64977 66029 66990 68861 69977 70037 71628 73336 76793 79488 
81891 81999 88003 90177 90263 90732 90749 94198 95638 97551 97733 97763 
97772 99145 99254 106033 106113 114278 115826 121308 121556 121789 123740 
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217313 


219134 
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229285 
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231394 
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236088 237605 


239525 
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240726 


240783 


242592 


242713 



242924 244065 248258 248482 248939 

5 -1 -2 (2): 952 26925 76782 84543 92035 94091 97702 105060 123639 
129074 145233 163463 206169 209432 217709 218420 222105 222396 223884 
229592 237971 242560 

5 -1 -1 (3): 97741 209474 

5 1-2 (4): 30887 47629 48023 48117 59592 73345 75496 118814 126994 
148536 149752 184749 
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5 2-2 (5): 3742 11237 48022 48138 48528 51225 75787 134216 147223 
150902 166364 173403 180497 190240 192315 221121 227919 

5 3 -3 (5): 41188 51919 52365 

5 3-2 (6): 28572 48012 48035 63365 69144 70259 77353 107672 118218 
147173 153032 173280 177352 180405 180485 180487 180488 180490 180495 
180496 180499 181592 195111 207078 208062 213634 217526 218932 220275 
220992 

6 -7 -1 (2): 9905 36976 45278 48016 48137 53029 56026 61235 64335 
65955 73093 94863 99656 128688 132181 134246 142870 145801 154396 158662 
167079 180492 181579 183958 185852 189036 199811 202232 205870 207658 
213096 214864 217403 224231 232466 

6 -7 5 (4): 21416 122051 180505 224163 228870 

5 4 -3 (6): 204141 

6 -6 -1 (1): 5684 60723 61210 77084 77140 78925 86835 92462 129754 
130529 135577 145896 158281 161636 164255 168790 169644 177138 177673 
192546 199074 216654 234556 

6 -4 -1 (1): 48001 97759 141744 157079 217839 

6 -3 -2 (1): 13797 19987 107393 122771 132844 204025 209370 237930 
244227 246410 

6 -3 -1 (2): 76800 

6 -1 -2 (3): 3426 11792 20861 91347 98826 142623 148819 170636 196499 200520 

6 1-3 (4): 57861 69553 76017 78776 116376 159801 168157 171677 
184224 202322 209447 209449 209451 209453 209462 209480 209484 219975 
232662 238386 240451 

6 1-2 (5): 16601 23644 30162 41621 68649 89802 92860 99716 121093 
132149 134702 139246 142066 150091 155432 163292 169468 188240 194745 
197727 197853 211543 215096 220219 228253 230149 241162 

6 2-3 (5): 4600 41401 51861 121579 159818 184675 209434 209467 
244168 

7 -7 -2 (2): 15290 15731 23485 47732 113231 170857 177020 177095 
181071 209450 212665 216400 236505 248433 
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7 -7 -1 (1): 218044 

7 -6 -2 (1): 9838 29563 97767 106034 127629 140618 178715 196816 
223929 248801 

7 -6 -1 (0): 33865 82020 

7 -5 -2 (0): 145459 

7 -4 -2 (1): 789 6473 9297 12276 30274 31230 51905 59439 62081 69966 71403 
78848 89247 93712 94018 96004 97368 112284 117217 126886 132654 143342 
150176 155000 178973 180648 181301 183437 184816 190375 207310 210646 
213978 215744 215955 216219 218726 225563 228100 232806 236561 237914 
238195 

7 -3 -2 (2): 16296 16431 37885 46792 47602 53099 61263 72777 88472 
93131 117093 119590 148796 162187 178526 183397 202167 219327 225991 
233905 239178 

7 -2 -3 (2): 1371 

7 -2 -2 (3): 10971 11148 16111 18757 20182 67775 70235 75374 95307 
95778 126511 153543 156313 166240 177296 180507 183679 194470 222451 
225756 228498 

7 -1 -2 (4): 20157 22519 27267 29032 33884 41473 48008 48027 48029 
50177 51198 53766 57387 62411 72909 79157 79676 88742 95068 99518 125118 
125567 130635 132085 132091 147946 154375 158303 159034 163148 170184 
174793 175160 175831 213528 214159 218114 240121 245749 

7 1 -3 (5): 51664 57935 65732 116950 190207 223715 

7 1 -2 (6): 48020 48031 68793 131749 174692 

7 2 -3 (6): 63304 66880 79094 145551 219356 237227 237964 240175 

7 3 -4 (6): 120291 

8 -7 -2 (1): 7282 12592 31488 36675 69219 72509 77518 99201 123051 
150021 157386 161417 166474 181601 182133 185382 186607 196547 204898 
206018 219402 232499 236599 244538 

8 -4 -3 (1): 10 468 51250 56857 113931 134378 173978 209448 209485 
8 -3 -3 (2): 105866 
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8 -3 -2 (3): 9791 21634 105959 106040 136430 181378 200353 225201 

8 -2 -3 (3): 151028 203738 

8 -1 -3 (4): 19698 41253 160914 185903 

8 -1 -2 (5): 181412 

8 1-4 (5): 221458 

8 1-3 (6): 4047 4553 7439 46087 96738 116854 125638 148867 152574 
219335 
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Appendix B 

Asteroids in pure two-body resonances with Jupiter 

1/1 (0): 588 617 624 659 884 911 1143 1172 1208 1404 1437 1583 1647 1749 
1867 1869 1870 1871 1873 2146 2148 2207 2223 2357 2363 2456 2674 2759 2797 
2893 2895 2920 3063 3240 3317 3391 3451 3540 3548 3564 3596 3708 3709 3793 
3794 4007 4035 4057 4060 4063 4068 4086 4138 4489 4501 4707 4708 4709 4715 
4722 4754 4791 4792 4805 4827 4828 4829 4832 4833 4834 4835 4836 4867 4902 
5012 5023 5025 5027 5028 5041 5119 5120 5123 5126 5130 5144 5209 5233 5244 
5254 5257 5258 5259 5264 5283 5284 5285 5436 5476 5511 5637 5638 5648 5652 
5907 6002 6443 6545 6997 6998 7119 7152 7214 7352 7543 7641 7815 8060 8125 
8241 8317 9023 9030 9142 9430 9431 9590 9694 9712 9713 9790 9799 9807 9817 
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249258 249465 249476 249481 249486 

2/1 (1): 1362 1921 1922 3688 3789 4177 5370 5528 8373 11097 11266 
11573 11665 13963 14871 16882 18888 22740 23577 26553 28368 28459 29524 
31249 31293 34901 35989 36140 37357 37528 37991 38406 38536 38584 38984 
39018 39309 41262 42202 45511 45796 46168 46204 47547 51070 51267 51537 
52700 55192 55222 55522 56326 57838 62679 62760 65297 65541 66174 68706 
68738 71661 71694 73396 73972 73995 76694 77869 77911 78154 78801 78814 
79282 79482 79500 82009 83904 83943 87362 87375 89175 89908 91182 91875 
95942 97406 97653 98988 99686 100100 102419 102915 104742 107690 113723 
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245143 245382 246128 

3/2 (1): 153 190 361 499 748 958 1038 1162 1180 1202 1212 1268 1269 
1345 1439 1512 1529 1578 1746 1748 1754 1877 1902 1911 1941 2067 2312 2483 
2624 2760 2959 3134 3202 3254 3290 3415 3514 3557 3561 3571 3577 3655 3694 
3843 3923 3990 4230 4255 4317 4446 4495 4757 5368 5439 5603 5661 5711 5928 
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6124 6237 6984 7027 7174 7284 8086 8130 8376 8550 8551 8743 8913 8915 9661 
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242451 


242656 


242682 


242851 


242996 


243267 


243311 


243438 


244397 


244438 


244478 


244488 


244878 


245007 


245093 


245216 


245578 


246178 


246496 


246497 


246578 


246726 


246786 


247405 


247507 


247553 


247563 


247854 


247879 


247893 


247900 


247920 


248004 


248430 


248945 


249017 


249060 


249182 


249246 


249266 


249275 


249284 


249302 



249324 249391 249416 

3/1 (2): 887 6318 48012 48016 48017 48022 48024 48032 48034 48035 
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48040 48138 48140 86924 92664 131702 144861 152664 180485 180487 180488 
180490 180492 180495 180496 180497 180498 180499 190119 207648 219527 
220909 

4/3 (1): 185290 186024 

5/2 (3): 26760 26817 48038 214027 230979 

5/1 (4): 209454 

7/3 (4): 99169 141224 182274 209430 209432 209433 209434 209435 
209439 209441 209443 209445 209452 209465 209466 209467 209470 

7/2 (5): 31010 48005 48006 48008 48010 48011 48014 48015 48020 48021 
48026 48027 48029 48030 48031 48033 48036 48037 48126 203386 

8/3 (5): 9671 30843 32859 48019 48023 48139 104777 120251 125664 
126492 129661 139227 150951 167506 168083 196568 210028 213384 226035 
237892 

9/5 (4): 209490 

9/4 (5): 33016 47371 52714 74469 101642 105466 121614 128854 134821 
150258 151436 163389 172417 173187 182682 186231 187691 187739 196684 
200129 201196 209429 209436 209437 209438 209442 209444 209446 209447 
209449 209456 209457 209458 209460 209463 209468 209471 209472 209473 
209475 209476 209479 209482 209486 209487 209488 209491 209495 209497 
230786 232988 237349 243711 245998 248481 

9/2 (7): 48001 175934 209455 209459 209474 209481 209492 209493 

10/3 (7): 48002 48004 48007 48009 48013 48018 48025 48028 48039 
48115 180486 180489 180491 180493 180494 

11/5 (6): 4439 7992 10032 18105 26484 34952 38277 48003 57445 64959 
99193 105840 106360 118403 120639 140934 146096 150539 162959 175111 
176464 179413 186287 197393 198525 206171 208965 209431 209440 209448 
209451 209453 209461 209462 209464 209477 209478 209480 209483 209484 
209485 209494 209498 211894 213992 240485 242608 247708 248231 

11/3 (8): 65720 98009 141845 161599 
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